THE ABLER KOSTANT SYMES SCHEME IN PHYSICS 



GABRIELA P. OVANDO 

Abstract. The purpose of this material is to review the Adler Kostant Symes 
scheme as a theory which can be developped succesfuUy in different contexts. It 
was useful to describe some mechanical systems, the so called generalized Toda, and 
now it was proved to be a tool for the study of the linear approach to the motion of 
n uncoupled harmonic oscillators. The complete integrability of these systems has an 
algebraic description. In the original theory this is related to ad-invariant functions, 
but new examples show that new conditions should be investigated. 



1. Introduction 

In this work we are interested in the use of Lie theory to understand some Hamilton- 
ian systems. For the study of completely integrable systems one needs to identify the 
following: i) the symplectic structure, which gives the system its Hamiltonian character, 
ii) first integrals or constants of motion, iii) action angle variables, and the computation 
of their evolution. Indeed this is a very difficult approach but it is possible for systems 
related to certain Lie groups. 

To this end there are several methods with a common idea: the realization of canon- 
ical equations on Lie algebras, or on orbits of a certain action or on symmetric spaces. 
These ideas appeared in the 70's and were developped, under other by several authors 
as Adler, Fomenko, Kostant, Mischenko, Olshanetsy, Perelomov, Trofimov, Symes, etc. 
(see for instance jAdT] |F-Mlj |F-M2j [F^ [KoT] [OP] [P] and [Sy] and their refer- 
ences) . 

In any method all of the above steps i) ii) iii) are reflected by algebraic circumstances. 
One needs a way for imbedding a certain Hamiltonian system into a Lie algebra, effective 
methods for constructing sets of involution and the proof of the full integrability of a 
wide family of functions in involution. 

In this chapter we are concerned with so called Adler-Kostant-Symes scheme, which 
brings together a mathematical framework with Lie theory but also consequences in the 
dynamics of the Hamiltonian system. This method was successful when studying some 
mechanical systems such as the rigid body or the generalized Toda lattice [Ad2 j |Ko2| 
|Sy| |R2| . In this setting the phase space of the Hamiltonian systems become coadjoint 
orbits represented on a Lie algebra and the functions in involutions are presented as 
ad-invariant functions. On the one hand for this kind of functions, the corresponding 
Hamiltonian systems become a Lax equation and on the other hand they are in involu- 
tion on the orbits. Whenever studying Poisson commuting conditions the ad-invariance 
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property can be replaced by a weaker one as in |R1] . In the framework of this theory 
what we need is a Lie algebra with an ad-invariant metric, a splitting of this Lie algebra 
into a direct sum as vector subspaces of two subalgebras and a given function. These 
algebraic tools were used with semisimple Lie algebras, where the Killing form is the 
natural candidate for the ad-invariant metric. 

However there are more Lie algebras admitting an ad-invariant metric. We shall 
examplify here how can be applied the theory for semisimple Lie algebras, and also for 
other ones, such as the solvable ones. For the general case one should see that any Lie 
algebra with an ad-invariant metric can be constructed by a double extension procedure, 
whose more simple application follows from M™. In this way one gets a solvable Lie 
algebra g, that results a semidirect extension of the 2n-|-l-dimensional Heisenberg Lie 
algebra f)„ and that can be endowed with an ad-invariant metric which is an extension 
of a non degenerate bilinear form on M^*^. But for other cases the resulting Lie algebras 
could be no semisimple and no solvable. 

In any case, the Lie algebra g splits naturally as a direct sum of vector spaces of 
two subalgebras. Looking at the coadjoint orbits of one of the Lie subalgebras, one 
gets Hamiltonian systems on these orbits and one can identify the original Hamiltonian 
system with one of these. In particular for the restriction of the quadratic corresponding 
to the ad-invariant metric we obtain a Hamiltonian system that becomes a Lax equation, 
whose solution can be computed with the Adjoint representation. 

As example we work out the Toda lattice and the linear equation of motion of n- 
uncoupled harmonic oscillators. The first one corresponds to a semisimple Lie algebra, 
and the second one is associated to a solvable one. Furthermore it is proved that the 
Hamiltonian for the last one is completely integrable on all maximal orbits. We notice 
that the functions in involution we are making use, are not ad-invariant and they do 
not satisfy the involution conditions of |R1] . 

The setting for the second example applies for quadratic hamiltonians. The Pois- 
son commutativity conditions we get for some polynomials can be read off in the Lie 
algebra sp{n) of derivations of the Heisenberg Lie algebra of dimension 2n-|-l In 
particular for the case of the motion of n-uncoupled harmonic oscillators we need a 
abelian subalgebra in the Lie algebra of isometries of the Heisenberg Lie group EI„, 
endowed with its canonical inner product. 

This is not surprising if we consider that symplectic automorphisms of the Heisen- 
berg Lie group produce symplectic symmetries of p-mechanical, quantum and classical 
dynamics for more general systems than the linear ones (see |Ki2| ). 

The appearence of the Heisenberg Lie algebra related to the motion of n-uncoupled 
harmonic oscillators is not so surprising. In fact, it is known that in quantum mechanics 
a good approach to the simple harmonic oscillator is through the Heisenberg Lie algebra. 
In dimension three this is the Lie algebra generated by the position operator Q = 
multiplication by x, the momentum operator P = and 1 with the only non trivial 

commutation relation 



[Q,^'] = 1 
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These operators evolve according to the Heisenberg equations 

dt dt 

An attempt to relate the classical mechanical system of the linear approximation 
of the motion of n-uncoupled harmonic oscillators was presented by the theory of p- 
mechanics, which makes use of the representation theory of the Heisenberg Lie group 
to show that both quantum and classical mechanics can be derived from the same 
source (see for instance | Kil] |Ki2| ). This theory contructs a more general setting that 
unifies both quantum and classical mechanics. The starting point for p-mechanics is 
the method of orbit of Kirillov |K1] |K2] . which says that the orbits of the coadjoint 
representation of the Heisenberg Lie group parametrise all unitary irreducible repre- 
sentations [F]. Thus non commutative representations are known to be connected with 
quantum mechanics. In the contrast commutative representations are related to clas- 
sical mechanics in the observation that the union of one dimensional representations 
naturally acts as the classical phase space in p-mechanics. In this theory the time evo- 
lution of both quantum and classical mechanics observables can be derived from the 
time evolution of p-observables, choosen as particular functions or distributions on the 
Heisenberg Lie group. 

These considerations allow to suppose that new applications of the Adler Kostant 
Symes scheme are possible and maybe it comes a new time to understand old me- 
chanical systems with new tools, which should be developped for these purposes. As 
an introduction to the topic one can find exceptional ideas in the books of Arnold, 
Abraham and Marsden, Ratiu and Marsden, etc. all of them classics in the literature 
concerning classical mechanics. 

The chapter is organised as follows: in the first part we present basic ideas concerning 
symplectic geometry. The second part is devoted to the Adler-Kostant-Symes scheme 
and the third part to the examples: on the one hand the Toda lattice with generalization 
in ( |Ko2| and |Sy| ), and on the other hand the systems corresponding to quadratic 
Hamiltonians on R^". 



2. Basic notions on symplectic manifolds 

In this section we present the basic elements to work with symplectic geometry. Some 
texts concerning this topic are [L-M| [CdSj . 

Let M denote a differentiable manifold. 

Definition 2.1. A 2- form on M, oo is called a symplectic form if da; = and is non 

degenerate for every p € M. 

The pair (M, w) is a symplectic manifold. 

It follows that the dimension of M must be even. 
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Example 2.2. Let M?^ be the usual euclidean space equipped with global coordinates 
xi, . . . , Xn, yi, • • • , Un- The 2-form given by 



a; 

1=1 

defines a symplectic form on R^"'. 

Note that if ( , ) denotes the canonical inner product on M^" and J the canonical 
complex structure 

'0 -I^ 
I 



J 



where / is the identity n x n matrix, then 

uj{X,Y) = {X,JY) 

Example 2.3. If (Mi,a;i) and (M2,W2) are symplectic manifolds, the direct product 
Ml X M2 is a symplectic manifold. 

Example 2.4. Coadjoint orbits. Let G denote a Lie group with Lie algebra q and let g* 
be the dual space of g. The coadjoint action of G on g* is defined as: 

g ■ ip = (f o Ad{g'^) geG,(p€g*. 

Notice that the orbit throught 93 is the set G • (p = {g ■ if : g G G} and the isotropy 
subgroup at ip is G^ = {g G G : ipoAd{g~^) = tp}; thus as usual one has G-ip = G/G^p. 
The action of G on g* induces an action of 3 on g* as 

X ■ if = —ip o ad(X) X e g,(p e g* 

that cames from the derivative 

— exp(tX) • = —93 o ad(X), 
dt It^o 

in other words 

X{^p) = -(^oad(X) 

is the infinitesimal generator induced byXGgat^jEg*. 

Any coadjoint orbit G ■ ip is a symplectic manifold with the 2-form 

uj^{X,Y) = -ip{[X,Y]), ipeg*,X,Yeg. 

called the Kirillov-Kostant-Souriau symplectic structure. 



Locally any symplectic manifold looks like the example (j2.2p above. This is a classical 
result of Darboux. 

Theorem 2.5. Darboux. Let {M,uj) be a symplectic manifold. For every p ^ M 
there exists coordinate system {U, {x^, . . . , x", y^, . . . , y")) such that p G U and = 
Y17=i dyi A dxi. 
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The symplectic 2- form uj bilds a isomorphism on TxM for x E M. In fact, since 
LVx ■ TxM X TxM is non degenerate, there is a hnear isomorphism Kx : TxM TxM 
defined by 

Kxiv){u) = uJx{u,v) = {-iyu;){u). 

Let : M — > M be a differentiable function, since its differential belongs to T*M, 
via the isomorphism above we get a vector field Xh given by 

(1) XH{x)=Kx{dHx), 

that is, Xh is the vector field on M satisfying 

v{H) = dH{v) =uj{v,Xh) 

and this is called the Hamiltonian vector field associated to the Hamiltonian function 
H. 

Example 2.6. For the standard symplectic structure on R^" the isomorphism Kx is 
given by KxV = —Jv, where J denotes the canonical complex structure on M^". 
Let H £ C°°(M^"), its associated Hamiltonian vector field is 

where VH is the gradient of with respect to the canonical inner product. 
Definition 2.7. The Hamiltonian system for a Hamiltonian H € C°°(M, w) is 

(2) x'{t) = XH{x{t)). 

Example 2.8. Let if be a smooth function on M^", the Hamiltonian equation is the 
classical one 

^' - 9l_ 
It' - 

ill ~ dxi 

where is actually Xi(t), that is it depends on t, for all i (and also for any yi). 
Example 2.9. On M^" a quadratic Hamiltonian is a smooth function as 

H{x) = —{Ax^x) for A symmetric linear map, 

which yields the Hamiltonian system 

(3) x' = J Ax 

In classical mechanic this system describes "small oscillations" , that is, it approximates 
the motion of a particle on or equivalently the motion of n uncoupled particles on 
M, near an equilibrium position. 

For instance the motion of n-uncoupled harmonic oscillators near an equilibrium 
position can be approximated with H a quadratic Hamiltonian as above by taking 
A = 1; therefore ([3]) becomes 

vm = -X^{t) 
where x{t) = (xi(t), . . . ,Xn{t),yi{t), . . . ,yn{t))- 
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In classical mechanics it is usual to name the coordinates as Xi the position coordi- 
nates and m as the velocity coordinates for every i = 1, . . . ,n. 

Definition 2.10. A diffeomorphism (p on a. symplectic manifold (M, w) is symplectic 
if (j)*u! = u). 

Recall that the Lie derivative on a smooth manifold M given as 

where X is a vector field on M with one parameter group tpt and T a tensor, satisfies 
the following identities 

Lx = ixd + dix 
Lxiy = ihxY + iyLx = i[x,Y] + iyLx 

For a proof see for instance |Wa] . 

Definition 2.11. A vector field X on a symplectic manifold (M, w) is symplectic if 
Lx^^ = 0. 

Proposition 2.12. A vector field X G x{^) is symplectic if and only if the one 
parameter subgroup tpt generated by X is symplectic. 

Proof. If ipt is symplectic, using the definition of Lx it is easy to see that Lx^^ = 0. 
Conversely assume Lx^^ = 0, then 

" ^1 ^^^^*'^ " Lxi'*sUJ = Tp*LxuJ = 
hence tp^uj is constant. But ipo = Id and so ip^uj = oo. □ 

Corollary 2.13. i) If lv is a symplectic form then Lx^^ = dix^- 

a) A vector field on (M, u) is symplectic if and only if ix^ is closed. 

Definition 2.14. A vector field X on a symplectic manifold (M, cj) is Hamiltonian if 
and only if —ix^ is exact. 

The vector field associated to a Hamiltonian function H defined in ([T]) is Hamiltonian. 

Notice that the fact of being X Hamiltonian says that there is H £ C°°{M) such 
that dH = —ixuj, therefore X is symplectic. On the other hand for any p £ M there 
always exists local solutions to dH = —ix^ for any X € x{^)- For global solutions 
we must ask extra conditions as below. 

Proposition 2.15. Let {M,uj) be a symplectic manifold such that i7-'^(M, M) = 0. 
Every symplectic vector field on M is Hamiltonian. 

A symplectic 2-form on a symplectic manifold M induces a Poisson bracket { , } 
on C°°M by: 

= ^p{Xf,Xg) = Xf{g) = -Xgif) for any f,g e C^M. 

Proposition 2.16. Let C°°{M) is a Lie algebra under the Poisson bracket defined 
above and f — > Xf is a Lie algebra anti-homomorphism of C°°{M) into x(M). 
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Proof. Since ■ T^M T*M is a linear isomorphism, the map / — > Xf is linear. 
Now we should prove that [Xf,Xg] = X^gj. Using the properties of Lx one gets 

LxfiXgUJ = i[Xj>,Xg]^^ + iXgLxfUJ. 

Since Lxft^ = 0, one gets 

LxfiXgUJ = i[Xf,Xg](^- 
The Lie derivative on 1-forms follows 

LxO = ixdO + dixO- 

Taking ixg^ = dg and applying above it holds 

Lxfixg^ = Lxfdg = ixfd^g + dixfdg = d{Xf{g)) = d{f,g}. 

Therefore 

HXf,Xg]^ = d{f,g} 

Thus the left side of the equality above i[Xf.Xg]^x coincides with —Kx{Xy^gj, and since 
Kx is an isomorphism = —X^jgj. □ 

Recall that a Poisson structure is a bracket { , } on a associative algebra A, such 
that 

• { , } is a Lie bracket on A and 

• f{g, h} = {fg, h} + {g, fh} for all f,g,he A. 

the last one is called the Leibnitz rule. In ISyl a such structure is called Hamiltonian. 



The space of smooth functions on a differentiable manifold is a associative Lie algebra, 
hence a natural space to be endowed with a Poisson structure. The Proposition we 
already proved says that whenever (M, w) is a symplectic manifold, C°°(M) has a 
Poisson structure induced by uj: the Poisson bracket { , } is a Lie bracket and the 
Leibnitz rule holds, since any vector field is a derivation on C°°(M). 

Example 2.17. On M^", the Poisson structure associated to the standard symplectic 
form is given by 

Example 2.18. Let g be a Lie algebra and q* its dual. As usual one identifies q* with 
its tangent space. Given a function F : g* ^ M, we define the gradient of F at a G g*, 
denoted by VF(a), as an element VF{a) £ g such that VF(a)) = dFa{P) for any 
/3 G g*, where ( , ) denotes the evaluation map. 
The Kirillov's Poisson bracket on g* is given by 

{f,h}{a) = {alpha, [V f{a),Vh{a)]). 

Proposition 2.19. If {f,g} = then g is constant on the integral curves of Xf. 

Proof. Assume x'{t) = Xf{x(t)) then 

|g(x(t)) = dgix'it)) = Xfig){xit)) = {/, <?}(x(t)) = 0. 

n 
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Thus g is called a constant of motion of the flow defined by Xf. Since { , } is skew 
symmetric, g is a constant of motion Xf if and only if / is a constant of motion of 
Xg. Constant of motion always exist, in fact / is a constant of motion of Xf. 

Definition 2.20. A function / on a 2n-dimensional Poisson manifold (M, { , }) is 
completely integrable if there exist n functions /i, . . . , /„ G C°°M such that: 

i) {/, fi] = 0, {fu /il = for all 1 < i, j < n, 

ii) The differentials dfi, . . . , dfn are linearly independent on a open set invariant under 
the flow of Xf. 

Two functions f,g : M ^ M such that {f,g} = are said to be in involution or 
Poisson commute. 

A subset C M is invariant under the flow of Aj if the solution x for the Hamil- 
tonian system ([2]) corresponding to the Hamiltonian / lies on A^ if x(0) € A^. 

Example 2.21. On M^" for H{x) = ^{x,x) the polynomials 

fi{x) = ^{pf + qf) i = l,...,n 

shows that H is completely integrable. In fact it is easy to check that {H,fi} = = 
{fijj} for ah i = 1, . . . ,n. 

Let F = (/i, . . . , fn), then F~^{c) is a torus which is invariant under the flow gen- 
erated by Xfj- Let {6i, . . . ,9n) denote the angle variable on the torus F~^(c). Then 
(/i, . . . , fn, Oi, . . . , On) is a local coordinate on R-^". With these coordinates, the Hamil- 
tonian equation becames 




and the coordinate functions satisfy 

{/i; fj} = %} = 0; %} = ^ij^ 

therefore the flow Xh is linear on F~^{c) for c € M". 

Moreover since the level sets {x G M^"" : H{x) = c} are compact we have action angle 
coordinates (see Liouville Theorem below). 

Generally m Poisson commuting functions /i , /2 , • • • , /m on a symplectic manifold 
(M, (j) give rise to an action of on M. Let {ipi)t be the one parameter subgroup 
generated by Xf.. Then 

{ti, . . . ,tm) ■ P = (^l)ti(V'2)t2 • • • 

defines a M*^ action on M. Since {fi, fj} = for all i,j, the set N = {x G M, : 
fi{x) = Ci} is invariant under the R^-action, for constants ci, . . . , Cm- If A^ is compact, 
the R^-action descends to a torus action on A^. When m = 1/2 dim M, one gets the 
Liouville theorem. 

Theorem 2.22 (Liouville). Let f he a completely integrable function on M, with 
dimM = 2n, and assume fi := /, /2,...,/n are commuting Hamiltonians which are 
linearly independent and let F = (/i, . . . , fn) ■ M — > R" he proper. Then F-^{c) IS 
invariant under the R" action and it descends to a torus T'^-action. Let 9i,...,6n 
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denote the angle coordinates on the invariant tori. Then {/«,/,} = {0i,6j} = and 
{fi,Oj} = Cij{F) for some functions Cij : M" — > M. In particular, the flow of Xf in 
coordinates (/i, ...,/„, ^i, ... , 9n) is linear. 

Coordinates as above, are called action-angle variables for the Hamiltonian system 
of/. 

3. SyMPLECTIC ACTIONS: THE AKS-SCHEME 

Let M denote a differ entiable manifold and let G be a Lie group. An differ entiable 
action of G on M is a differentiable map rj : G x M ^ M, ry : (g, m) rj{g, m) := g - m 
such that 

i) e ■ m = m for all m € M and 

ii) {gh) ■ m = g ■ {h ■ m) for all m € M, g,h ^ G. 

Notice that if rj is an action, the applications rjg : M ^ M given by r]g{m) = g ■ 
m are diffeomorphisms of M. In fact, rjg are differentiable for any g and they are 
diffeomorphisms since the inverse of any rjg is 77^-1 (see ii) above). Therefore an action 
of a Lie group on M induces a representation of G on Diff(M) the diffeomorphisms 
of M, given hy g ^ rjg. 

Example 3.1. Let GL{n,M) denote the Lie group of non singular transformations of 
M". This acts on M" as evaluation: A ■ v = v io'c A ^ GL{n,'R) and G M". It is easy 
to verify that this is in fact an action. 

Example 3.2. Let H he a Lie subgroup of a Lie group G, then H acts on G by conju- 
gation, H X G —>■ G, {h,x) = h^^xh, for any h ^ H , x ^ G. If is a normal subgroup, 
one can consider the action of G on by conjugation. 

Example 3.3. Let G be a Lie group with Lie algebra g, then G acts on g by the Adjoint 
action, G x g ^ g, (^r, X) = Kd{g)X, for any g G G, X G g. Recall that A-d{g) = dl{g)e 
where Ig denotes the conjugation by g (see the previous example). It is easy to see 
that Igh = IgO Ij^ for all g,h £ G, hence the map G GL{q) is a representation of G, 
called the Adjoint representation. This has a correlative at the Lie algebra level, the 
adjoint representation: g x g ^ g given by X -Y = [X, Y] for all X,Y £ g. 

Recall that in (j2.4p we defined the coadjoint action of a Lie group G on the space 
g*, the dual of the Lie algebra 

g ■ Lf = (f o Ad{g'^) 5 G G, G g* 

and also we gave the corresponding action of g on g* by 

X ■ (f = —ip o ad(X) X G g, G g*. 

The orbit of an action of a Lie group G on a set M is 

G ■ m = {g ■ m : g £ G} 

and the isotropy or stabilizer group of the action at the point m is the closed subgroup 
of G given by 

Gm = {5 G G such that g ■ m = m}. 
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It is known that the orbit at m is diffeomorphic to the quotient space of G and the 
isotropy group, G • m G/Gm (see |Wa] for instance). Thus any curve at the orbit 
G ■ m through m is ^{t) = exptX • m and this generates the infinitesimal vector X at 
r^(G-m) by 

Xim) = — exritX ■ m. 

Hence the tangent space of a G-orbit at m is 

Tm{G-m) = {X, 

being g the Lie algebra of G. 

Assume M and are two differentiable maps on which a given Lie group G acts. 
A map F : M ^ A^ is called equivariant if F{g ■ m) = g ■ F{m) for all m € M, g £ G. 
The condition is also expressed as F intertwines the two G-actions. 

Definition 3.4. Let {M, to) be a symplectic manifold. An action 77 of a Lie group G 
on M is called symplectic if the diffeomorphisms rjg are symplectic maps for any g G, 
that is r]*ui = lo. 

The coadjoint orbits are examples of symplectic manifolds. Recall that they are 
endowed with the 2-form given by: 

uj(,{X,Y) = -P{[X,Y]), PeG-fi 

which is symplectic. In fact, it is closed since for Xi,X2,X3 G g one has 

uj{[Xi,X2],Xs) = -yp{[[Xi,X2],Xs]]), 

hence 

dcu(Ai, A2, A3) = -(^([[Ai, A2],X3]]) - (^([[A2,X3],Xi]]) - (^([[Aa, Ai], A2]]) = 

where the last equality holds after Jacobi for [•,•]. 

The 2-form uj is non degenerate on a orbit: let ip £ q* and let A G g such that 
u;{X,Y) = for all Y £ g. 

Then -ip{[X,Y]) = for all Y £ q, says that X ■ ip = implying that A £ L{G^). 
In fact exptAT £ G^ if and only if expiA •</? = </? for t near 0. Thus taking derivative 
at t = we have AT • 93 = 0, and this is the set corresponding to the Lie algebra of G^p. 
Since the tangent space of the orbit at ip is r^(G • ip) = q/L{G^), one gets X = 0. 

Definition 3.5. An ad-invariant metric on g is a bilinear map ( , ) : g x g — > M, which 
is a non-degenerate symmetric and such that ad(A) is skew symmetric for any X £ q, 
that is 

([A, Y],Z) + {Y, [A, Z]) = for ah X,Y,Z £ g. 

This ad-invariant metric gives rise to a bi-invariant pseudo Riemannian metric on a 
connected Lie group G with Lie algebra g; bi-invariant means that the maps Ad(5f) are 
isometries for all g £ G, that is 



{Adig)Y,Adig)Z) = {Y,Z) 



for all Y,Z £ Q,g £ G, 
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and conversely any bi-invariant pseudo Riemannian metric on G induces an ad-invariant 
metric on its Lie algebra, just by taking derivative of the last equality at t = with 
g = exptX. 

Examples of Lie algebras with ad-invariant metrics are: 

a) semisimple Lie algebras with the Killing form; 

b) semidirect products g t<coad9* with the canonical neutral metric 

An ad- invariant metric ( , ) induces a diffeomorphism between the adjoint orbit G X 
and the coadjoint orbit G ■ ix where ix{Y) = {X,Y). In fact 

g ■ £x{Y) = {X, Ad{g-^)Y) = {Ad{g)X, Y) for all X,Ye0,geG, 

implying that the map i : X ^ ix is equivariant. Thus the adjoint orbits become 
symplectic manifolds with the 2-form: 

iOx{Y,Z) = {X,[Y,Z]) for X,y,ZG0. 

Wc shall consider these ideas to construct Hamiltonian systems on orbits that are 
included on Lie algebras. 

Recall that given a metric ( , ) on g the gradient of a function / : g ^ R at the vector 
X G g is defined by 

(6) {Vf{X),Y)=dfx{Y) yeg. 

Suppose g+, g_ are Lie subalgebras of the Lie algebra g such that 

g = g+©g_ 

as a direct sum of linear subspaces, that is (g,g+,g_) is a product structure on g. The 
Lie algebra g also splits as 

g = gi©g^, 

and 

g^ is isomorphic as vector spaces to g^. 

This follows from the isomorphism ^ : g ^ g*. In fact, let X G gij: maps to ix- 
Since ix{Y) = for all Y e g+, the image of £{q^ belongs to g^l, and the isomorphism 
follows from dimensions. 

Let G- denote a subgroup of G with Lie algebra g_. Then the coadjoint action of 
G- on g!. induces an action of G_ on q^: for g^ G G_, X G gij:, F G g_ one has: 

g--ix{Y) = {X,Ad{9-^)Y) = {Ad{g)X,Y) = 7T^j_{Ad{g.)X),Y), 

where tt^x denotes the projection of g on giji; therefore the action is given as 

5--X = 7rgx(Ad(5-)X), 

and i : ^ g'^ is equivariant. 

The infinitesimal generator corresponding to y_ G g- is 

Y^{X) = -^ expty_-X = 7r X G g|. 
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The orbit G_ • Y becomes a symplectic manifold with the symplectic structure given 
by 

ojx{U-,V-) = {X,[U-,V-]) for U-,V- £Q-,X eG--Y 

which is induced from the Kostant-Kirillov-Souriau symplectic form on the coadjoint 
orbits in q*_. 

Consider a smooth function / : g ^ M and restrict it to an orbit G- ■ X := A4 <Z 
Q^. Then the Hamiltonian vector field of the restriction H = f\j^ is the infinitesimal 
generator corresponding to — V/_ , that is 

(7) XH(y) = -vrgx([v/_(y),y]) 

where Z± denotes the projection of Z € with respect to the decomposition g = 



g+ © g_. In fact for Y G g:}:, V- G g_ we have 



uy{v.,Xh) = d/7y(y_) = (v/(y),7rgx([F_,y])) = (v/_(y),[F_,y]) 
= (y,[v/_(y),y_])=^y(v/:(y),F_). 

Since to is non degenerate, one gets ([7]). 

Therefore the Hamiltonian equation for x : M ^ g follows 

(8) x'(t) = -7rgx([V/_(x),x]). 

In particular if / is ad-invariant then = [Vf{Y),Y] = [V/_(y),y] + [Vf+{Y),Y]. 
Since the metric is ad-invariant [g+,g:j:] C g:}:, in fact 

([0+,0i],0+) = (si, [0+,0+]) = 0. 

Hence the equation ([8]) takes the form 

(9) x'{t) = [VU{x),x\ = [x,Vf.{x)l 

that is, ([8|) becomes a Lax equation, that is, it can be written as x' = [P{x), x\. 

If we assume now that the multiplication map x G- — > G, {g+,g-) — > g+g-, is 
a diffeomorphism, then the initial value problem 



(10) 



f = [V/+(x),x] 

x(0) = Xq 



can be solved by factorization. In fact if exptV/(xo) = g-\-{t)g^{t), then x{t) = 
Ad((7_|_(t))xo is the solution of ([10]). 

Remark. If the multiplication map G+ x G- — > G is a bijection onto an open subset 
of G, then equation ([8]) has a local solution in an interval {—e,e) for some e > 0. 

The theory we already exposed shows the application of Lie theory to the study of 
ODE's as in equation ([9]). Even when it is possible to give the solution, one need more 
information. This can be obtained from involution conditions. They help in some sense 
to control the solutions. 

A first step in the construction of action angle variables is to search for functions 
which Poisson commute. The Adler-Kostant-Symes Theorem shows a way to get func- 
tions in involution on the orbits A4. We shall formulate it in its classical Lie algebra 
setting. 
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Theorem 3.6 (Adler-Kostant-Symes). Let q he a Lie algebra with an ad-invariant 
metric ( , ). Assume 0-, 0+ are Lie subalgebras such that g = g_ g+ as direct sum of 
vector subspaces. Then any pair of ad-invariant functions on g Poisson commute on 
g+ (resp. on q-^). 

Sometimes the ad-invariant condition is too strong, so the following version of the 
previous Theorem given by Ratiu [Rl] asks for a weaker condition. 

Theorem 3.7. Let q he a Lie algebra carrying an ad-invariant metric ( , ). Assume it 
admits a splitting into a direct sum as vector spaces g = g+ ©g_, where g+ is an ideal 
and g_ is a Lie subalgebra. If f,h are smooth Poisson commuting functions on g, then 
the restrictions of f and h to g^ are in involution in the Poisson structure o/g^. 

Remark. This theorem was used in [R2] to prove the involution of the Manakov integrals 
for the free n-dimensional rigid body motion. 

4. Applications of the Adler-Kostant-Symes-scheme to classical 

mechanics 

In this section we show the explicit use of the theory above in some Lie groups and 
Lie algebras. The first example is done with semisimple Lie algebras, and it is known 
as the Toda Lattice. 

4.1. The Toda lattice. The Toda lattice is the mechanical system which describes 
the motion of n particles on a line with an exponential restoring force, that is the 
Hamiltonian function on M^n is 

^ n n—1 
i=l i=l 

The phase space is M^" which is a symplectic manifold with its canonical symplectic 
structure. It follows that the Hamiltonian equation is 



(11) ^'^ ~ 

and with e^°~^^ = = e^"~^"+i. Flaschka considered a change of coordinates (called 
Falschka transform) as follows 

0:M2n^]^2n^ ^(a;, y) = (a, 6) 

where 

flfc = -\yk I <k <n, 

hk = le^^'k-^k+i l<k<n-l 

bn = ^6 2 



Therefore the equation (jlip yields 



b'k = bk{ak+i - fflfc) I <k <n 
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with an+1 = = 6o- Notice that Y^iUi = 0. Assume X^a;^ = = and let 

V = {{x, y)/ Ylii = = Ylii yi}j then the system above becomes 

(12) < = 2(^f.-&ti) l<k<n-l 

b'k = h(.ak+i -ak) 1 < /c < n - 1, 

Consider g the semisimple Lie algebra of traceless real matrices sl{n,M. equipped 
with the ad-invariant metric {x,y) = tr{x,y) for all x,y (z s/(n,M). 

Let 0+ = so(n) the Lie subalgebra of skew symmetric real matrices and 0_ the Lie 
algebra of upper triangular matrices of trace zero. 

Then gij: is the space of real symmetric matrices in sZ(n,M) and g-^) is the space of 
strictly upper triangular matrices in sl{n,M). 

The coadjoint orbit ^A = G--xo for xq = "^IZi + is the set of tri-diagonal 
real symmetric matrices 

n n—1 

aiEi^i + ^ bi{Ei^i+i + Ei+i^i) ^ ai = 0, 6i > Vi. 

i=l i=l i 

where Eij denotes the matrix with a 1 at the place i,j and in the others components. 

Let / ■ sl{n,R) ^ R be the function given by f{X) = ^{X,X) = ^tr{XX). It is 
easy to see that the gradient of / at X is X, and hence applying the theory of the 
previous section we get ([9]) 

x' = [x+, x] 

for x+ C 5+ a curve in Writing the last system in terms of coordinates {a,b) we 
get the system ([T2]). 

A generalization of this system can be read in [Syj . where also aplications of the 
theory to other differential equations are explained. 

4.2. The motion of n uncoupled Harmonic oscillators. Recall that the motion 
of n-uncoupled harmonic oscillators near an equilibrium position can be approximated 
with H the quadratic Hamiltonian as ^{x, x) where ( , ) is the canonical inner product 
in M^". Let a; the canonical symplectic structure, the corresponding Hamitonian system 
follows 

where x{t) = . . . , x„(t), . . . , 

The associated Poisson structure on M^" is given as follows 

df dg df dg 



^^'^^ ^■^'^^ (V/, Q^.Qy, dyidxi' 

for smooth functions /, g on R^". Thus with respect to the Lie bracket { , } the subspace 
over R generated by the functions H = ^ + yi )^ the coordinates Xi, yi, and 1 

form a solvable Lie algebra of dimension 2n+2, which is a semidirect extension of the 
Heisenberg Lie algebra spanned by the functions Xi,yi, 1 i=l, . . .,n. In fact they obey 
the following non trivial rules 

{xi,yj} = 5ij {H, Xi} = -yi {H, yi} = Xi. 
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In order to simplify notations let us rename these elements identifying Xn+i with 
H, Xi with Xi, Yi with yi and Xo with the constant function 1 

1 ^ Xo 

Xi <— Xi 
Vi ^ Yi 

H Xn+i 

and set g denotes the Lie algebra generated by these vectors with the Lie bracket 
[• , •] derived from the Poisson structure. This Lie algebra is known as a oscillator Lie 
algebra. 

The Lie algebra Q splits into a vector space direct sum g = g_|_ ©g_, where q± denote 
the Lie subalgebras 

(15) 0- = span{XQ, Xi, Yj}ij=i^,„n, 0+ = ^Xn+i- 

Notice that Q- is isomorphic to the 2n+l-dimensional Heisenberg Lie algebra we 
denote 

The quadratic form on g which for X = xo(X)Xo + + Vii^Wi) + 

Xn+i{X)Xn+i is given by 

/(^) = \ ^i^i + Vi) + ^OXn+l 
i 

induces an ad- invariant metric on g denoted by ( , ) . It is easy to show that the gradient 
of / at a point X is 

V/(X) = X. 

The restriction of the quadratic form to B := span{Xi,Yj} i, j=l, . . ., n, coincides 
with the canonical one ( , ) on R^" ~ t>. 

The metric induces a decomposition of the Lie algebra g into a vector subspace direct 
sum of and g^ where 

= span{Xo} = MXn+i © span{Xi, y,-}ij=i,...,n, 

and it also induces linear isomorphisms gj_ — g^. Let G denote a Lie group with Lie 
algebra g and G± C G is a Lie subgroup whose Lie algebra is g±. Hence the Lie 
subgroup G- acts on by the "coadjoint" representation; which in terms of U- G g_ 
and F G gij; is given by 

(16) adlrV = Xn+i{V)Y,i{yi{U)Xi-Xi{U)Yi) 

It is not difficult to sec that the orbits are 2n-dimensional if Xn+iiV) and further- 
more V and W belong to the same orbit if and only if Xn+i{V) = Xn+i{W), hence the 
orbits are parametrized by the x„+i-coordinate; so we denote them by Mx^.^.^- They 
are topologically like M^". In fact Mx„+^ = G_ • V ~ EI„/Z(EI„), where EI„ denotes 
the Heisenberg Lie group with center Z(]H[„). 

Equipp these coadjoint orbits with the canonical symplectic structure, that is for 
G g_ take 

n 

u;y{U-,V.) = {Y, [U., V.]) = Xn+i{Y) Y,ixi{U-)yi{V.) - XiiV.)yi{U.)). 

1=1 
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Indeed on the orbit M.i the coordinates Xi,yj, i,j = l,...n, are the canonical sym- 
plectic coordinates and one can identify this orbit with M^" in a natural way. This 
says that the identification is a symplectomorphism between M^" with the canonical 
symplectic structure and the orbit with the Kirillov-Kostant-Souriau symplectic form. 

Consider H, the restriction to a orbit Aix^+i of the function /. Since / is ad-invariant 
the Hamiltonian system of = /|^ reduces to 



(17) 



dx 
dt 

x{0) 



[Xn+lXn+l,Xx) + Xn+lXn+l] 

n.0 



where a;° = + x^^^Xq and = + ViYi)- 

For Xn+i = a^n+i = 1 t^i^ system is that one we get on 

The trajectories x{t) with coordinates Xi{t), yj{t), x^_^i are parametrized circles of 
angular velocity x^^^, for all i,j, that is 



Xi{t) 

yjit) 

Xn+l{t) 



a;° cos(4 , li) + y° sin(x°+it) 
-x^ sin(4+ii) + cos{xl^^t) 



X, 



n+l 



This solution coincides with that computed in the previous section, when we considered 
systems on coadjoint orbits. In fact it can be written as 

x{t) = Ad{exp tx'^^iXn+i)x^ , 

and one verifies that the flow at the point X° G is 

A*(X°) = EJ(^?cos(xO+it) + i/Osin(4+,t))X, + (-xOsin(xO+ii)+ 

(18) 

Vi cos(x°+it))Fi] + 
By taking L and M the following matrices: 



M 








Xn+1 

















-Xn+1 






























Xn+1 


















-Xn+1 


























Xn+1 
















-Xn+1 




























v 

















0/ 
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L = 






Xn+l 














Xl\ 


Xn+l 

















Vi 











Xn+l 








X2 



















y2 












a^n.+i 





Xn 






-^2/2 


\X2 .. 


-Xn+l 








yn 

































0/ 



we get L' = [M, L] = ML — LM, the Lax pair equation. 



5. Quadratic Hamiltonians and coadjoint orbits 

In this section we shall prove that Hamiltonian systems corresponding to quadratic 
Hamiltonians in M^" of the form H{x) = \{Ax, x) where A is a non singular symmetric 
map, can be described using the scheme of Adler-Kostant-Symes on a solvable Lie 
algebra. 

Let us consider the linear system of one degree of freedom on M^" with Hamiltonian 
given by: 

H{x) = ^(Ax,x) 

where x = {qi, . . . ,qn,pi, ■ ■ ■ ,Pn) is a vector in M^" written in a symplectic basis and A 
is a non singular symmetric linear operator with respect to the canonical inner product 
( , ) . This yields the following Hamiltonian equation 

dSD x' = J Ax, with J = 

and being Id the identity. The phase space for this system is M^". We shall construct 
a solvable Lie algebra that admits an ad-invariant metric on which the system ([3]) can 
be realized as a Hamiltonian system on coadjoint orbits. Moreover it can be written as 
a Lax pair equation. 

Let b denote the non degenerate bilinear form on M?^ = span{Xi,Yj}'^j^^ given by 
b{X,Y) = {AX,Y). In our terms, b defines a metric on M?^ but it is not necessary 
definite. Note that the linear J A is non singular and skew symmetric with respect to 
b, where J is the canonical complex structure on ti ~ M^" as above: 

b{JAX, Y) = {AJAX, Y) = {J AX, AY) = -{AX, JAY) = -b{X, J A). 

Let Q denote the Lie algebra g which as vector space is the diract sum 

Q = mxo e D e mx„+i 

where = M^" and with the Lie bracket given by the non trivial relations 

(19) [U, V] = b{JAU, V)Xo [Xn+l, U] = JAU for all [/ G 0. 

Thus in this way one defines a structure of a solvable Lie algebra on q. Note that 
^ = Id is the particular case we considered in the previous subsection. 
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This Lie algebra q can be equipped with the ad-invariant metric defined by 

(20) {xlXo + U^+xi+iXn+i,xlXo + U^+xl+^Xn+i) = 6([/i,C/2) + (xJ4+i+x24^i). 

Thus if (, )p denotes the restriction of the metric of g to D = spanjXj, i^}jj=i^...^„, 
then clearly ( , ) is a generalization of the non degenerate symmetric bilinear map b of 
M^". Moreover q admits a orthogonal splitting 

£1 = span{Xo, Xn+i} t). 

Denote by g± the Lie subalgebras 

0+ = RXn+i, 0- = MXq © span{Xi, , Yi}. 

They induce the splitting of into a vector space direct sum = 0+ ©0_, which by the 
ad-invariant metric gives the following linear decomposition = 0:}: ® 0:*: , direct sum 
as vector spaces, for 

01^ = RXq 0^ = span{Xi, yOi=i,...,n © RX^+i- 

Note that 0_ is an ideal of isomorphic to the 2n-|-l-dimensional Heisenberg Lie algebra 

Let G denote a Lie group with Lie algebra 0, set G_ C G the Lie subgroup with Lie 
subalgebra 0_. As we already explained G_ acts on by the coadjoint action 

5--X = 7rgx(Ad(5-)X) g-eG., X e g^, 

where tt^j. is the projection of on 0:}:, which in infinitesimal terms gives the following 
action of 0_ on 

(21) axi*uV:=U-V = Xn+i{V)JAX^{U) for C/ G 0_, F G 0:^. 

being X^{U) the projection of U onto B with respect to the orthogonal splitting = 
span{Xo, Xn+i} © 0. 

The orbits are 2n-dimensional if Xn+i (V) ^ and furthermore V and W belong to the 
same orbit if and only if Xn+i{V) = Xn+iiW), and therefore one parametrizes the orbits 
by the a;„+i-coordinate and one enotes them by M-x^+i ■ The orbits are topologically 
like M^" since they arc diffeomorphic to the quotient ]HI„/Z(1HI„), if Z{Mn) = MXq is the 
center of the Heisenberg subgroup. 

Endow the orbits with the canonical symplectic structure of the coadjoint orbits, 
that is for X E Q^, U-,V- e Q- set 

ux{U-,V.) = {X, [U.,V.]) = Xn+i{X)b{JAU^,V^). 
Consider / : — M the ad-invariant function given by 

f{X) = ^{X,X). 

The gradient of the function / at a point X is the so called position vector 

VfiX)=X. 
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Since / is ad-invariant the Hamiltonian system of H = f\ , the restriction of / to 
the orbit M.x„+i , given by Q becomes 



dx 

(22) 



-j^ - [V/+(x),x] = [Xn+lXn+l,Xo + Xn+lXn+l] = Xn+l-JAx^ 

x(0) = 



where X^ G g^. 

Thus this Hamiltonian system written as a Lax pair equation is equivalent to Q for 
j^i — 'f- The solution X{t) for the initial condition X^ G g:^ can be computed 
via the Adjoint map on G, that is, 

X{t) = Kd{exp 

The previous explanations prove the following result. 

Theorem 5.1 ( |Q2j ). Let H{X) = ^{AX^X) he a quadratic Hamiltonian on M^" with 
corresponding Hamiltonian system Then H can he extended to a quadratic function 
f on a solvable Lie algehra g containing the Heisenberg Lie algebra as a proper ideal. 
The function f induces a Hamiltonian system on coadjoint orbits of the Heisenberg 
Lie group, that can be written as a Lax pair equation and which is equivalent to 
Moreover the trajectories on M?"^ for the initial condition can be computed with help 
of the Adjoint map on g. Explicitely they are the curves x{t) = exp^'^^V^, where exp 
denotes the usual exponential map of matrices. 

If we take L, M G M(2n + 2, M) as 

(Xn+lJA z\ /xn+lJA 

L = I i\z'^ 
00/ \0 0^ 

where z'^ = (xi, X2, • • • , Xn, yi,y2, ■ ■ ■ , y-n) then the Hamiltonian equation can be written 
in the following way 

L' = [M,L]. 

Example 5.2 (The motion of n-uncoupled inverse pendula). As example of the previous 
construction consider the linear approximation of the motion of n uncoupled inverse 
pendula. This corresponds to the Hamiltonian H{x) = ^{Ax,x) with 

Ld 
-Ld 



This yields the Hamiltonian system x' = J Ax, which in coordinates takes the form 

dx 

X. 



(23) 4 

dt 

As we said the phase space is M^". In the setting of the AKS scheme we can construct 
coadjoint orbits M of the Heisenberg Lie group, that are included in a solvable Lie 
algebra g with Lie bracket (llOh and ad- invariant metric (I20p . The Hamiltonian system 



20 



GABRIELA P. OVANDO 



for the restriction to the orbits of the ad-invariant function on q, f{X) 
be written as 

[Xn+lXn+l,Xx, + Xn+lXn+l] 



h{X,X), can 



(24) 



dx 
dt 

x{0) 



where X^ = + Vi^) + ^n+i^n+i- The Hamiltonian system above on the 

coadjoint orbit Adi written in coordinates is clearly equivalent to (j23p . 

The trajectories on q^, x = Ylii^ii't)-^i + yi{'t)^i) + Xn+iXn+i are parametrized by 

Xi{t) = x^cosh(x°^_it) + sinh(x°^_it) 
yi{t) = x°sinh(x°+it) + y°cosh(x°+it) 

Xn+l{t) = X^^i 




The flow at the point X G is 
(25) 



A\X^) = ^.[(xOcosh(:EO 



n+it) - y°sinh(j;0+^t)Xi+ 



+(a;° sinh(x°+it) + y° cosh(x°+it)yi] + xS^iX, 



The system ()24p is a Lax pair equation L' = [M, L] = ML — LM , and has a matricial 
representation by choosing L and M the following matrices in M(2n + 2,M): 



M 







Xn+1 





















Xn+1 


































Xn+1 






















Xn+1 






























Xn+1 




















Xn+1 







































V 
















0/ 









Xn+1 














Xl 


\ 


Xn+1 

















yi 





















X2 










Xn+1 











y2 



















Xn 










^X2 


Xn+1 








Vn 




Ixi 


■ ■ ~2^" 










V 


























/ 



Now we shall investigate involution conditions on the coadjoint orbits of the Heisen- 
berg Lie group for the restrictions of the quadratic functions f{X) = ^(X, X), where 
( , ) denotes the ad- invariant metric on the solvable Lie algebra g. 

Let gi-iQj be two quadratics on M^" that are realted to the symmetric maps Ai,Aj : 
t) ^ t) respectively, that is 



1 



g,{X) = -{AiX,X) 



gj{X) = -iAjX,X). 
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Consider quadratic functions on the solvable Lie algebra g, which are extensions of 
gi,gj to RXo ® MX^+i, for instance as 

1 1 
gi{X) = -{AiX^,X„) + xoXn+i 9j{X) = -{AjX^,X„) + xox^+i. 

For the following results these extensions are not unique. For instance extending them 
trivially we get the same conclusions. 

Let Hi,Hj denote the restrictions oi gi,gj to the orbits J^xn+\ and let X G C 
g^. The symplectic structure on the orbits induces a Poisson bracket which for the 
functions Hi,Hj follows: 

{Hi,H,}{X) = {X, [Vg,_{X),Vg,_{X)]) 

By computing one can see that the gradients of gi and gj are 

Wgi{X) = A-^AiX„ + xoXo + Xn+iXn+i Vgj{X) = A~^AjX„ + xoXo + Xn+iXn+i. 
Thus we are ready to prove the following result. 

Theorem 5.3 ( |Q2| ). The functions Hi,Hj are in involution on the orbits Mx^+i ^/ 
and only if 

(26) [JA„JAj]=0 
where J is the canonical complex structure on M^". 

Proof. Let X E M-x^+i C 0+. For the functions Hi, Hj the Poisson bracket on the orbit 
follows: 

{H,,Hj}{X) = {X,[AiX„AjX,]) = {xn+^[X^+^,A-^AX,],A'^AjX,) 
= Xn+i{JAiX„,A-^AjX„) = Xn+i{JAiX„,AjX„) 

Therefore {Hi, Hj}(X) = if and only if (AjJAiX^j,X„) = which is equivalent to 
AjJAi = AiJAj, if and only if JAjJAi = JAiJAj, that is [JAi, JAj] =0. □ 

The natural question is what is the meanning of (j26p ? 

Fix ( , )' the inner product on f)„ defined so that the vectors Xi,Yj,XQ are orthonor- 
mal for all i,j = l,. . ., n. The metric is an extension of the canonical one on M^". The 
Lie bracket on [)„ = RXq © where M^" ~ o = spanjXj, l^}jj=i^...^„ is expressed as 

{[X,Y],xoXoy = xo{JX,Yy with J as in ^ 

and note that (, )|^^^ = (, ). A derivation D of f)„ acting trivially on the center must 
satisfy [DU,V] = —[U,DV] for all U,V G v. Equivalently in terms of (, )', we have 
that a map D in [)„ is a derivation acting trivially on the center of [)„ if and only if the 
restriction of D to (denoted also D) satisfies 

{JDU, V) = -{JU, DV) for all U,V 

where we replaced ( , )(, by ( , ) since they coincide on ~ M^". Denote by t) the set of 
derivations on f)„ acting trivially on the center of 

Theorem 5.4. There is a bijection between the set of derivations of f)„ acting trivially 
on the center and the set 5o{n) of symmetric linear maps on M^". This correspondence 
is given by D d ^ JD € so(n), where J is the complex structure as in 
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Corollary 5.5. If there exists an n- dimensional ahelian suhalgehra on z{JA)x), where 

z{JA)^ = {D ed such that [D, J A] = 0} 

then the Hamiltonian function H restriction of the function f{X) = ^{AX, X) is com- 
pletely integrable on the orbits Mx^+i Z^?^ ^n+i 7^ 0. 

Proof. The previous theorem says that the restrictions to the orbit of the func- 

tions gi, gj are in involution if their corresponding derivations commute in t). In particu- 
lar for gi and /, we have that H and Hi Poisson commute on the orbit if and only if JAi 
belongs to the centralizer of J A in d, z{JA)^. Since the complete integrability requires 
of n linearly independent functions, this can be done with a basis of an n-dimensional 
abelian subalgebra of z{JA)i,, finishing the proof. □ 

A linear map t is a derivation of P)„ acting trivially on the center }{\^)n) if and only 
if Jt + t*J = 0, if and only if t € sp(n). The derivations of nilpotent Lie algebras of 
H-type were computed in ([Sii]). 

In the case of the motion of n-uncoupled harmonic oscillators, we can see that the 
corresponding derivation is an element of a Cartan subalgebra of sp(n). 
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